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ABSTRACT 


The  variation  of  the  reactivity  of  a  pile  as  a  fimction 
of  width  of  a  transverse  gap  is  obtained.  The  method  involves 
first  finding  the  boundary  condition  satisfied  by  the  flux  at 
the  gap  face.  This,  in  principle,  provides  enough  information 
for  a  complete  solution  of  the  pile  equations.  A  method  for 
calcixlating  the  reactivity  change  is  presented.  The  calculated 
reactivity  is  compared  with  experiment  and  a  brief  discussion 
of  the  validity  of  the  approximations  is  given. 
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THE  EFFECT  OF  GAPS  ON  PILE  REACTIYITT 
S .  Tamor 

I.  Introduction 

The  problem  of  a  transverse  gap  in  a  cylindrical  reactor  has  been 
treated  in  some  detail  by  Goldberger,  Goldberger,  and  Wilkins  several 
years  ago^^ .  They  develop  a  method  for  calculating  the  flux  depression 
at  the  gap  vhich  is  somewhat  restricted  in  its  validity  because  of  an 
approximation  which  is  invalid  for  small  gaps.  The  method  used  here  is 
exactly  that  of  GGW,  but  the  approximation  to  large  gaps  is  not  made. 

No  serious  complications  is  introduced  by  the  generalization.  Once  the 
boundary  condition  at  the  gap  is  obtained,  it  is  possible*  to  replace 
the  gap  by  an  equivalent  thin  absorber;  i.e.,  an  absorber  which  produces 
the  same  fliix  depression.  It  is  then  a  simple  matter  to  calculate  the 
reachivity  change  by  means  of  perturbation  theory  (assuming,  of  course, 
that  the  adjoint  f l\ix  for  the  reactor  in  question  is  known) . 

At  the  risk  of  repeating  previous  work,  for  the  sake  of  complete¬ 
ness  the  principles  of  the  GGW  method  will  be  presented  in  detail.  We 
will  restrict  ourselves  to  reactors  in  the  shape  of  a  right  (not  neces¬ 
sarily  circular)  cylinder.  As  pointed  out  by  GGW  the  general  method  is 
applicable  to  other  geometries,  the  only  requirement  being  that  the  flux 
be  separable  into  parts  normal  and  transverse  to  the  gap. 

II.  The  Method 

Consider  a  bare  cylindrical  reactor  of  length  which  has  been  cut 
along  a  transverse  plane,  and  the  two  parts  separated  a  distance  h 
(fig.  1).  For  convenience,  we  let  the  gap  lie  in  the  mid-plane  of  the 


1) 


Goldberger,  Goldberger,  and  Wilkins,  CP  hereafter  referred  to 

as  GGW . 
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reactor  although  the  generalization  to  an  off-center  gap  is  trivial. 
Ttie  flux  in  the  one—  velocity  diffusion  problem  is  separable  into 
transverse  and  longitudinal  parts,  so  that  within  the  moderator  the 
longitudinal  part  is  a  cosine  distribution  eind  the  transverse  part 
satisfies  the  two  dimensional  wave  equation  'V  0  +  70  =  0  with 

appropriate  boundary  conditions. 


—  2  - > 

Figure  1 


We  assume  the  reactor  to  be  sufficiently  large  that  diffusion 
theory  is  applicable  and  gives  correctly  the  flux  emerging  from  the 
gap  face.  This  restriction  implies  that  for  reactors  with  several 
gaps  they  must  not  be  too  close  together.  Following  GGW,  we  now 
inquire  as  to  how  many  neutrons  leave  an  element  of  area  at  p '  on 
the  +  face  and  arrive  at  p  on  the  -  face. 

If-TL  is  a  unit  vector  directed  from  p '  on  the  +  face  to  p  on 
the  -  face,  the  rate  at  which  neutrons  leave  p '  headed  into  dlX  is 


-'(0.(P')  +  A.  V0.(^')'ix)cos  4  dil 

ifjt  + 


so  that  the  total  number  of  neutrons  arriving  at  p  is  given  by 


n  - 


4« 


0^  (p')  +  X  V  0^-  (p')-n 


cos  y  d/l 


-L 

4« 

where  A '  is  the  element  of  area  containing  p ' . 
to  the  total  incoming  current  at  p  which  is 


0;^(P')  +  X  V0^(P')^ 


cos 


<l'dA' 


(1) 


This  is  to  be  equated 


0.(^  X  c)0_(p) 

“V" 

Making  use  of  the  symmetry  of  the  system  about  the  mid-plane, 


.  )4.(?)  and  |!i'(p)  s  ?i,'(p)  E  3i'p) 

*  *  i  z .  '  ■  TV 


we  have 


^  ^  t-U' 


1 

4« 


0  +  X 


cos 


Ip  dA 


R 


(2) 


(3) 


At  this  point  it  is  desirable  to  make  one  more  approximation. 

Let  us  assume  the  gap  to  be  small  compared  to  the  transverse  dimensions 
of  the  reactor  so  that  there  is  no  appreciable  communication  between  the 
outside  and  the  middle  of  the  reactor.  This  is  the  region  of  maximum 
importance,  and  the  incoming  flux  will  not  be  seriously  falsified  by 
continuing  the  integration  in  (3)  over  the  infinite  plane.  We  will  see 
that  this  device  leads  to  a  very  great  simplification  in  that  the  spatial 
distribution  of  the  entering  flux  reproduces  that  of  the  flux  emerging. 
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To  evaluate  the  intergrals  it  is  convenient  to  view  figure  1  as 
projected  onto  the  gap  face  (fig.  2).  The  coordinate  system  is  chosen 
arbitrarily  and  u  is  directed  from  p  to  p ' 


Figure  2 


Here  the  GGW  calculation  neglects  the  last  term  in  the  numerator 
of  the  integrand.  It  is  easily  seen  that  this  is  equivalent  physically 
to  ignoring  the  azimuthal  asymmetry  of  the  emerging  flux.  Since  the 
flux  is  directed  preferentially  toward  regions  of  low  density,  the 
azimuthal  distribution  favors  those  neutrons  heading  toward  the  edge 
of  the  reactor,  thereby  contributing  a  net  streaming  out  of  the  gap. 

We  will  see  that  this  gives  most  of  the  neutron  loss  for  sufficiently 
small  gaps. 

Since 

a  u  ^ 


u 
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and  cos  =  h/R;  eq.  (4)  becomes 


q(p)  = 


4it 


(h  +u  )  I  Vu^+h2  '/u^+h2 


^0(p) 


h  dA'  (5) 


But  the  flux  is  separable  so  that  0'(p>)  is  equal  to  ^(p')  times  a 
function  of  z  which  is  constant  over  the  gap  face.  Therefore  the 


K(u)0(p')dA'  where  K(u)  is 


integral  may  be  written  in  the  form 
explicitly  a  function  of  u  alone. 

Choosing  the  point  p  as  a  new  origin,  the  Fourier  Bessel  trans¬ 
form  of  0(p')  may  be  written 


0(p’)  =  z 

t 


<is  a.g^(p)  J^(su)cos(tX  + 


Substituting  into  +  7^0  =  o  and  making 
property,  we  find 


use  of  the  orthogonality 


(s  -7  .  0 


SO  that 


0(p’)  =  ^  (p)J^(7u)cos(tX+ 

0(p')<iX 


(6) 


Row 

K(u)  0(p’)dA'  = 

P 

K(u)u  du 

n  ^ 

*  c 

- 

2«a^Q(p)cos  K(u)Jq(7u)u  du 

(7) 


But  we  now  observe  that  if  p 


(6)  becomes 


0(p)  =  a^Q(p)cos 

K(u)0(p')(3A'  = 


U 


so  that 


2it0(p) 


K(u)Jq(7u)u  du 


This  proves  our  assertion  that  the  incoming  flux  is  proportional  to 
the  emerging  fluxj  and  further  that  the  constant  of  proportionality 
depends  upon  h,  y,  but  is  independent  of  the  geometry  of  the  system. 
Returning  to  the  original  integral,  (5)  becomes 


q(p)  = 


0(p) 


u  du-  h^7>^ 


(  2  ,  2,,2 

(u  +h  ) 


t  2  ,2^5/2  "  ^ 

(u  +h  Y' 


0(7u) 


u  du. 


These  integrals  are  all  Hankel  transforms  and  are 


2) 


Jo(7u) 

,  2  ,2-.  2 
(u  +h  ) 


u  du  = 


'0 


Jo(7u) 


7K^(  7h) 
2h 


u  du  = 


7^^^K3/2(7h) 

2^/\3/2p(5/2) 


7e 


•7h 


3h" 


^  Jj_(7u)  2  7^^^Ki/2(7h)  7e"’'‘ 

22/3tjl/2  p  (5/2)  3h 


2) 


See  Watson,  Bessel  Functions  p.^3^' 
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where  the  K  fiinctions  are  the  usual  Bessel  fvmctions  of  the  second 
kind  with  imaginary  argment.  We  now  have 

q(p)  =  ^  ■ 

^  f|.(s)2%2ii  e-’'" 


2  1  -  7kK^(rh)  +  I  (7X)(rh)e"’'^ 

^  1  +  (l+7h)e"^^ 


Note  that  as  h-»0  0  0  and  as  h  c»  ^  2X  ’  usual  diffusion 

theory  boundary  condition. 

The  result  of  GGW  is  identical  with  (9)  "but  for  the  omission  of 
the  last  term  in  the  n^^merator.  For  small  gaps,  the  last  term  is 
linear  in  h,  while  1  -  2Kj^(z)  starts  off  as  z  log  z.  Thus  for  very 
sma3_L  h  the  last  term  dominates.  For  reasonable  dimensions  of  the 
reactor,  the  two  parts  become  comparable  for  h  55J\/3.vhich  for  graphite 
is  about  1  cm.  For  gaps  much  larger  than  this  the  last  term  may  be 
neglected. 

For  small  h  eq.  (9)  can  be  checked  by  an  elementary  method.  Let 
h  be  small  enough  that  the  flux  and  its  angular  distribution  are 
essentially  the  same  in  the  gap  as  in  the  neighboring  moderator.  The 
total  fltix  leakage  out  of  the  gap  is  then 


I*-" 
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vhere  ds  is  an  element  of  arc  along  the,  gap  edge^  and  ^  is  the  out¬ 
ward  gradient  of  0  at  the  edge. 

The  rate  at  which  neutrons  are  fed  into  the  gap  is 


-  2 


V  0-djJ 


where  dA  is  an  element  of  area  on  the  +  face.  For  yh.  1 


(11) 


0'  y\ 
f  ^  2 


so  that  (11)  becomes 


S 


0dA. 


But  0  satisfies  V  0  +  7  0  =  0  giving 


-f 


0dA. 


Integrating  by  parts,  if  we  let  s,  be  a  unit  vector  along  ds,  we  have 

2  -a  ^ 

-  \h  (  Vx  s,  )0.ds.  Since  0  is  constant  along  the  gap  edge,  this 
^  “  2  ■*'r  dd 

is  simply  -  ^  ^  which  is  precisely  the  result  (10). 


III.  Calculation  of  the  Reactivity  Change 

So  far  we  have  considered  the  effect  of  the  gap  on  neutrons 
diffusing  with  constant  velocity.  While  this  is  an  unrealistic 
picture  of  a  reactor,  our  result  does  provide  enough  information  to 
obtain  the  6k  of  a  real  system.  For  a  tr\aly  thermal  reactor  the 
following  somewhat  crude  method  may  be  used. 
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We  start  with  the  pile  equation 
V  _  k  (B^) 

®  i+l¥ 

where  k  is  the  infinite  medium  multiplication  constant,  p  is  the 
resonance  escape  probability,  is  the  Fourier  transform  of 
the  infinite  meditim  slowing-down  kernel,  L  is  the  diffusion  length, 
and  B  the  geometrical  buckling. 

For  h  =  0  the  boundary  condition  at  the  gap  in  0 '  =  0 .  Knowledge 

0' 

of  ^  thus  determines  the  increase  in  buckling  of  the  system.  If  B 
is  small,  we  have 


5k 


eff 


eff 


6(B^) 


when  r  is  the  mean  square  distance  traveled  in  slowing  down  to 
thermal. 

Fitting  a  cosine  distribution  to  the  new  boundary  condition,  we 

find 


5(B^) 


k  0' 

It 


so  that 


k 


4 

J 


^  0' 


(12) 


3) 


See  A.  M.  Weinberg  and  L.  C.  Noderer,  Theory  of  Neutron  Chain  Reactions, 
Vol.  II,  Part  I,  Chapter  V,  (CF  51-5-98). 


This  method  assumes  that  all  the  fissions  are  caused  by  thermal 
neutrons  and  breaks  down  if  an  appreciable  fraction  of  the  fissions 
are  fast. 

It  is,  however,  possible  to  obtain  the  6k  for  more  general  reactors 

provided  that  the  flux  and  its  adjoint  for  the  unperturbed  system  are 

known.  This  can  be  done  by  noting  that  the  gap  may  be  replaced  by  an 

4) 

equivalent  thin  absorber  ' .  The  presence  of  the  gap  causes  a  net  flow 

2 

of  neutrons  into  the  gap  from  the  two  faces  of  ^  ' .  These  neutrons 

are  effectively  absorbed"  by  the  gap.  If  ^  is  small,  we  replace  the 
gap  by  an  absorber  whose  cross-section  is  Z(z)  =  9  6(z-z^)  where  z^ 
defines  the  plane  at  which  the  reactor  is  opened.  This  will  absorb  the 
correct  number  of  neutrons  if 


0Z(z)dz 


u 


or 


9 


006(z-z^)dz 


U 


0(Zo)0 


(13) 


The  effective  absorption  cross-section  is  seen  to  be  a  fimction  of 

energy  through  its  dependence  on  \. 

In  the  spirit  of  the  raultigroup  method  we  consider  the  slowing- 

down  process  as  a  succession  of  one-velocity  diffusion  proccesses  each 

providing  a  source  for  the  next  lower  one  in  energy.  At  each  energy 

the  result  (9)  is  applicable.  If  we  have  solved  the  slowing-down 

6k 

problem  for  the  unperturbed  reactor,  ^  may  be  obtained  from  perturbation 
theory.  We  have 


q''^(5^,u)dr'  du 


(14) 


4) 


This  method  was  suggested  to  the  author  by  D.  K.  Holmes. 
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where  q(r,u)  is  the  slowing-down  density  at  lethargy  u  and  coordinate 
r,  Z(z)  is  defined  by  (13) ^  and  the  other  symbols  have  their  usual 
meaning. 

For  a  bare  reactor  q(r,u)  is  separable  into  q(u)*Q(r)  where 
Q(r)  is  self-adjoint  and  (l4)  becomes 


6k 

“IE 


Q^(r)  b{z-z^)6T 


q(u)  q'^Cu)  du  , 


and  if  Q  is  properly  normalized 


8k 

k 


cos^6Zo)| 


q(u) 


e(u) 


q^(u)  du 


(15) 


where  6  is  the  longitudinal  buckling. 

It  is  easily  shown  that  (15)  is  valid  also  for  the  case  of  the 
off-center  gap.  The  result  implies  that  the  6k  produced  by  a  gap  of 
given  width  is  proportional  to  the  importance  of  the  plane  along  which 
it  is  cut. 

Substituting  into  (15)  from  (9)  and  (13)  we  have 


k 


where 


a  = 


2/^  N  2 

cos  (e  z^)  j 


1  -  (7h)K3_(7h) 
1  +  (l+7h)e"’'^ 


al. 


+  pig 


(17) 


P  = 


72h  e-’'^ 


^  1  +  (H-7h)e“’'^ 


In  = 


<l(u) 


q'‘’(u) 


^2  = 


The  GGW  result  is  obtained  by  setting  p  =  0. 

By  this  method  6k  is  calculable  in  a  straightforward  way  involving 
only  two  Integrations.  In  many  cases,  \  is  independent  of  u  over  a 
wide  range,  and  may  be  considered  constant  so  that  actually  only  one 
integral  need  by  evaluated. 

Validity  of  Approximations 

The  approximations  made  in  obtaining  the  boundary  conditions  at 
the  gap  are  a)  use  of  diffusion  theory  for  the  outgoing  flux  at  the 
gap  face,  and  b)  the  continuation  of  the  surface  integrals  to  infinity. 

A  qualitative  justification  of  the  diffusion  approximation  is 
simply  this.  The  diffusion  approximation  gets  poorer  as  the  singular 
distribution  of  the  neutron  flux  deviates  more  and  more  from  isotropy. 
The  worst  possible  case  from  this  point  of  view  is  when  h  becomes 
infinite.  However  this  leads  to  the  elementary  diffusion  boundary 
condition  which,  for  large  systems,  is  an  excellent  approximation  to 
the  correct  one.  The  error  arising  from  this  source  is  then  probably 
no  worse  than  that  from  using  the  diffusion  approximation  in  the  pile 
or  multigroup  equations. 

To  justify  the  continuation  to  infinity,  one  observes  that  the 

2  2  2  2 

integrals  involved  go  something  like  h  /(u  +  h  )  times  a  relatively 
slowly  varying  function  of  u.  The  first  factor  is  very  sharply  peaked 
in  the  neighborhood  of  u  =  h  and  falls  off  as  a  high  power  of  u.  If 
the  distance  from  the  edge  of  the  reactor  is  large  compared  to  h,  the 
controbution  to  the  integral  from  large  u  is  small.  For  points  close 
to  the  edge  the  approximation  breaks  down,  but  these  regions  have  low 
importance  and  csm  be  ignored.  The  error  was  estimated  by  numerical 
integration  at  the  center  of  a  circular  cylinder  of  radius  =  lOh 
(a  relatively  large  gap).  In  this  case  the  leakage  is  overestimated 
by  about  1.5  7®*  For  smaller  gaps  it  is  easy  to  see  that  the  accuracy 
improves  very  rapidly. 
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The  replacement  of  the  gap  by  an  equivalent  absorber  and  the 
use  of  perturbation  theory  require  that  the  flxix  depression  at  the 
gap  be  small. 

V.  Comparison  With  Experiment 

The  methods  developed  above  were  used  to  calculate  tVie  reactivity 

change  as  a  function  of  gap  width  for  the  CA-4  graphite  moderated 

criticeO.  assembly^ ^ .  This  reactor  was  a  rectangular  parallelepiped 

130  X  112  X  112  cm  with  a  gap  perpendicular  to  the  long  dimension  and 

6k 

slightly  off  center.  Measurements  of  for  gap  widths  up  to  about 
1  cm  were  made,  the  maximum  opening  determined  by  the  available  control. 

It  was  unfortunately  extremely  difficult  to  obtain  accurate  measurements 
of  the  table  separations  for  very  small  gaps,  hence  the  large  experimental 
uncertainty. 

Figure  3  compares  the  experimental  points  with  theory,  curve  1 
representing  the  GGW  result  obtained  from  eq.  (17)  but  with  3=0,  while 
curve  2  is  calculated  from  the  complete  expression.  Eq.  (12)  does  not 
apply  to  this  system  since  a  large  fraction  of  the  fissions  are  caused 
by  fast  neutrons.  For  this  case  eq.  (12)  gives  about  one-third  of  the 
total  effect. 

In  view  of  the  rather  large  uncertainty  in  the  data  the  check 
between  theory  and  experiment  should  be  taken  as  only  qualitative.  At 
this  point  it  cannot  be  said  whether  the  apparent  divergence  between 
theory  and  experiment  for  large  h  is  real. 

The  comparison  is,  however,  sufficiently  encouraging  that  for 
reactor  in  which  diffusion  theory  is  applicable  one  may  expect  to  predict 
the  effect  of  a  small  gap  with  some  degree  of  confidence. 

5) 

For  further  details  of  the  experiment  see  Report  Y-881. 
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FIGURE  3.  REACTIVITY  LOSS  vs  TABLE  SEPARATION 

I  (f  =  ^  OF  7.3  X  lO""^ 

^  k 
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